带非线性支撑的转子有限元模型求解方法

用数值方法研究了非线性支撑的柔性转子系统的动学行为,提出了一种将有限元与非线性支撑结合的模型和求解方法。利用有限元法(FEM)构建转轴和转盘部分的模型,通过矩阵进行组合;利用离散元方法对包含滚动轴承和挤压油膜阻尼器(SFD)的支撑部分进行建模,此部分包含4个单元,分别为轴承内圈、外圈、SFD内圈和支撑鼠笼。有限元部分和离散元部分通过轴端节点相连,仿真过程中轴端位移传递给非线性支撑部分,支撑部分通过位移计算得到的非线性力反过来作用于有限元转子轴端部分。为了耦合求解有限元转子和非线性支撑组成的数学模型,提出了一种综合的迭代求解方法,克服传统的有限元求解方法对轴端隐性非线性支撑的求解局限性。由于转轴部分采用了Timoshenko梁单元建模,对比与简单转子模型,可以考虑陀螺力矩和轴的柔性特征,更能体现非线性支撑对振动真实影响。在建立的20个轴单元的有限元转子模型中,非线性响应更多体现在靠近非线性支撑的节点1和节点21处,响应频谱中靠近轴端的节点能体现出滚动轴承的2倍和3倍变柔振动频率。     

水下超声速射流对上浮水雷受力特性影响研究

固体火箭发动机具有功率密度大、推力大等优势,常被用于上浮水雷的推进器。然而发动机在水下工作时燃气射流使流场压力发生剧烈脉动,进而影响发动机的推力性能及上浮水雷的受力特性。基于VOF多相流模型和理想气体模型,建立了上浮水雷在不同工作状态(欠膨胀、完全膨胀、过膨胀)下的数值模型,研究了水下燃气喷射流对上浮水雷的受力特性影响。结果显示,欠膨胀工况时,发动机推力平缓,大小为12.2KN,上浮水雷受力未出现负值;完全膨胀时,射流发生颈缩、胀鼓现象位置距离雷体较远,导致发动机推力及水雷壳体受力振荡不剧烈;过膨胀工况时,射流发生颈缩、胀鼓现象距离喷管较近,发动机推力发生剧烈脉动,产生21.37%的振荡幅度,胀鼓现象发生时,流场压力显著降低使得上浮水雷后体受力减小,壳体阻力增大,上浮水雷最大产生27KN的负推力。     

基于音圈电机的环形桁架结构主动悬吊方法

环形桁架结构是卫星可展开天线的基本结构。文章提出一种主动悬吊方法,可在地面状态下卸载环形桁架结构的重力。该方法基于音圈电机可以同时输出安培力和大位移,通过求解悬吊系统的静态方程计算悬吊点的张力即所需安培力;然后,建立主动悬吊试验系统,令音圈电机输出等于所需安培力,即实现悬吊结构对重力的卸载。经试验验证,悬吊后的结构前两阶频率与原结构的几乎相等,且模态一致,说明该悬吊方法有效。     

推力室逆置点火瞬时真空舱流场数值与试验研究

针对高空模拟舱内逆向布置的推力室点火瞬时过程中舱内压力的变化规律进行了理论和数值仿真分析,推力室高速气流的动能是推力室短程点火过程中舱内压力大幅上升的主要原因.为保证推力室点火时舱内的真空度,需要在推力室出口配置导流装置将燃气导出真空舱.对单台和4台并联推力室进行点火试验,试验表明：采用燃气导流能够较好保持舱内真空度,是整机多推力室高空模拟试验的新思路.     

含孔边裂纹各向异性有限板的二维问题

应用边界元法，并结合一种特殊的复变函数基本解，研究了含孔边裂纹的各向异性有限板的应力强度因子修正系数的计算。这种特解边界元法避免了在裂纹附近划分单元，可以获得较好的数值计算精度。通过算例分析，结果表明本方法精度高，运算量小，具有一定的工程应用价值。     

反推状态下大涵道比涡扇发动机气动稳定性预测与评估

为了预测与评估反推状态下,反推气流再吸入对大涵道比涡扇发动机气动稳定性的影响,采用反推气流扰流流场三维CFD数值模拟、发动机整机稳定性计算分析以及反推状态下发动机进气畸变台架试验相结合的方法,开展了反推气流对大涵道比涡扇发动机气动稳定性影响的研究。通过三维CFD数值模拟手段,捕获了反推状态下发动机进口流场的畸变程度。在此基础上,通过采用稳定性计算程序预测了发动机的气动稳定性,并进一步通过发动机台架试验,验证了预测结果。CFD计算结果表明,随着相对来流马赫数的减小,反推气流被发动机重新吸入的可能性不断增大,当相对来流马赫数减小到0.05时,外侧发动机进口的流场畸变情况变得最为严重。进气畸变情况下的整机稳定性计算分析以及发动机台架试验结果表明,在所考核的目标状态,若只存在因反推气流再吸入引起的进口流场畸变,是不会导致发动机失稳的。     

我国可重复使用液体火箭发动机发展的思考

重复使用是降低航天发射成本的重要途径之一,是液体火箭发动机未来发展的重要方向。本文分析了可重复使用液体发动机的发展趋势,针对可重复使用运载器对发动机功能的需求,探讨了动力系统方案;对比了液氧煤油和液氧甲烷等推进剂组合和不同循环方式,认为几种发动机方案均可满足重复使用运载器的需求;研究了重复使用发动机的关键技术,提出应重点研究可重复使用液体火箭发动机高温组件热结构疲劳寿命评估及延寿技术、运动组件摩擦磨损技术、结构动载荷控制与评估技术、快速检测评估与维修维护技术、健康监控与故障诊断技术、二次或多次起动技术与大范围推力调节技术等。     

基于终端角度约束的二阶滑模制导律设计

针对空地导弹具有终端角度约束条件的制导律设计问题,提出了一种在有限时间内稳定的新型二阶滑模制导律。首先,在弹目相对运动学模型基础上,将终端弹道倾角约束转化为终端视线(LOS)角度约束,作为制导系统的终端控制目标。其次,通过选取一种新型二阶滑模面,结合螺旋控制算法的思想,设计了一种二阶滑模变结构制导律,来抑制系统中的不确定性因素,从而满足零化视线角速率和制导系统的终端角度约束条件的要求。采用一种新的Lyapunov函数,基于Lyapunov稳定性理论,严格证明了制导系统在有限时间内的稳定性。最后,对空地导弹制导系统进行数字仿真,通过和一阶传统滑模制导律以及基于超螺旋算法的二阶滑模制导律进行对比分析,验证了所设计的制导律在保证制导精度的同时,更能在有限时间内提高终端约束角度的精度,并且避免了超螺旋算法中参数选取较多的问题。     

翼型测力天平的研制

为研究翼型结冰状态下的载荷变化规律，对翼型进行测力实验。研制一台三分量翼型测力天平。天平设计的突出难点为设计载荷极不匹配，阻力测量载荷很小。设计时，针对翼型载荷的特点，采用片梁组合方式测量阻力，四柱梁测量法向力及俯仰力矩，并采用有限元方法对结构进行了验证计算。主要介绍天平的设计、校准及应用情况。     

小波变换数值均匀化方法求解振荡系数的方程

将小波数值均匀化方法用于求解二维快速振荡系数的抛物型方程.小波数值均匀化方法基于小波的多分辨分析,在细尺度空间上得到原方程的半离散方程,然后利用小波映射得到粗尺度空间的数值均匀化方程,即可在粗尺度上求解原方程,有效地减小了计算量.计算结果表明:小波数值均匀化方法与精细剖分的有限差分法相比较,既大大地节省了计算时间又获得了较好的精度.